timesGPU= [24 51 79 150];
timesGPU2= [31 41 59 101];
timesREF= [103 500 920 1997];
points = [10000 50000 100000 200000];
speedups = timesREF./timesGPU;
speedups2 = timesREF./timesGPU2;

figure(1);
hold on;
xlabel('Points of integration');
ylabel('Time (s)');
set(gca,'yLim',[0.0 2100.0]);
set(gca,'xLim',[0.0 210000.0]);
set(gca,'YTick', timesREF);
set(gca,'YTickLabel', timesREF);
set(gca,'XTick', points);
set(gca,'XTickLabel',points);

scatter(points,timesGPU,'r','^','filled');
scatter(points,timesGPU2,'g','v','filled');
scatter(points,timesREF,'b','o','filled');
plot(points,timesGPU,'r');
plot(points,timesGPU2,'g');    
plot(points,timesREF,'b');
legend('C2050','K40c','Original','Location','Best')

figure(2);
hold on;
xlabel('Points of integration');
ylabel('Speedup');
set(gca,'yLim',[0.0 25.0]);
set(gca,'xLim',[0.0 210000.0]);
set(gca,'YTick', [3     4    10    12     13    16    20]);
set(gca,'YTickLabel', [3     4    10    12     13    16    20]);
set(gca,'XTick', points);
set(gca,'XTickLabel',points);
scatter(points,speedups,'r','^','filled');
scatter(points,speedups2,'g','v','filled');
plot(points,speedups,'r');
plot(points,speedups2,'g');
legend('C2050','K40c','Location','Best')